Time evolution of tetragonal-orthorhombic ferroelastics 
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We study numerically the time evolution of two-dimensional (2D) domain patterns in proper 
tetragonal-orthorhombic (T-O) ferroelastics. Our results, found by solving equations of motion 
derived from classical elasticity theory, disagree with those found by other methods. We study first 
the growth of the 2D nucleus resulting from homogeneous nucleation events. The later shape of the 
nucleus is largely independent of how it was nucleated. In soft systems, the nucleus forms a flower- 
like pattern. In stiff systems, which seem to be more realistic, it forms an X shape with twinned 
arms in the 110 and 1 10 directions. Second, we study the relaxation that follows completion of 
the phase transition; at these times, the T phase has disappeared and both O variants are present, 
segregated into domains separated by domain walls. We observe a variety of coarsening mechanisms, 
most of them counterintuitive. 



PACS numbers: 81.30.Kf, 68.35.-p, 62.20.Dc, 61 
I. INTRODUCTION 
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Ferroelasticara are crystalline solids that undergo a 
shape-changing phase transition, usually first-order, to 
a state of lower symmetry with decreasing temperature 
T. A prominent example is the-tetragonal-orthorhombic 
(T-O) ferroelastic YBa 2 Cu3O7.0t!l At the transition tem- 
perature T c , the unit cell of the parent (high-T) phase 
distorts spontaneously in one of several equivalent direc- 
tions. Each of these degenerate distortions corresponds 
to a differently oriented variant of the product (low-T) 
phase. Below T c , all variants are usually present, sepa- 
rated from each other by domain walls with preferred ori- 
entations. Domain patterns in ferroelastics differ greatly 
from those in ferromagnets, gainsaying.-lhe analogy re- 
sponsible for the very name ferroelastiuH and confound- 
ing intuition based on conventional order-parameter sys- 
tems. The difference from these other systems is that the 
strains are not independent order parameters but rather 
are linked by compatibility relations. 

The theory of proper ferroelastics (where the strain 
is the primary order parameter) extends the classical 
theory of elastic continua by adding higher-order terms 
in the strains and also derivatives of the strains. This 
strain-only theory was first usedo in one dimension ClD); 
its first important result was a remarkable solutions for 
the twin wall of cubic-tetragonal (C-T) materials. It 
has since been used to study various aspects of (a) T-O 
materials, ul3 (bJ-the-lD problem,Ej (c) cubic-tetragonal 
(C-T) materials, Eiil3 and (d) hexagonal-orthorhombic 
(H-O) and related materials.113 Although the strain-only 
theory applies strictly only to proper ferroelastics, it has 
nevertheless succeeded in explaining_a wide variety of do- 
main patterns also in improper T-OE3 and H-OO materi- 
als. A much larger literature (examples are Refs. |l^|2^) 
includes order parameters in addition to the strains or 
applies more phenomenological approaches. 

The basic formalism for describing the time evolution, 
of proper ferroelastics, though known for a century,^ 
has been used only infrequently, to study 1D,Ej C-Tcf 



and H- systems. Fundamentally different dynami- 

cal schemes were used in the strain-only theories of Refs. 
[7],p| pl|^5pl| , often only as a tool to find static structures. 

The following presentsjthe first application of the clas- 
sical equations of motionEj to the dynamics of proper T- 
O ferroelastics. The study .ajas motivated in part by the 
electron-microscopy resultsou available for YBa2Cu307; 
this is an improper material (the orthorhombic distor- 
tion is a secondary effect of the oxygen-ordering) , but the 
success of the static strain-only theoryO for YBa2Cu307 
warrants an extension to the dynamics. Computational 
resources allow us to consider only 2D structures, with 
possible application to thin films, particularly to the pat- 
terns of Refs. UH 

The paper is organised as follows. Section II gives the 
expression for the T-O strain energy in 2D and then finds 
the equations of motion. We distinguish between soft and 
stiff systems according to the energy cost for wall direc- 
tions off optimal. The strain-only theory predicts soft- 
ening with decreasing T, perhaps with observable conse- 
quences. Section III applies this formalism to investigate 
the growth of O nuclei from the supercooled T phase. 
We find that the developed nuclei are largely indepen- 
dent of the nucleation mechanism; in both soft and stiff 
systems rtlaey differ markedly from the nuclei in other 
theories. OE3 In soft systems, nuclei are flower-like; they 
simply expand without generating much additional struc- 
ture. In stiff systems, they form an X shape with twinned 
arms in the major growth directions (110 and 110); ad- 
ditional structure forms near the centre and propagates 
outward along the arms. Section IV examines the coars- 
ening mechanisms that follow completion of the phase 
transition, including domain-wall merges, formation and 
disappearance of island domains, rank formation of rib- 
bon tips and their coordinated retraction, and tip split- 
ting (in stiff systems). Section V provides a summary 
and proposals for further investigations. 
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II. EQUATIONS OF MOTION 

A. Expansion of the strain energy 

The energy of proper ferroelastics is expressed solely in 
terms of the strains. These are combinations of deriva- 
tives of the displacement u(x) of a material point from 
its position x in the high-T symmetric phase. We dis- 
cuss only structures uniform in the tetragonal fourfold 
(£3) direction. We define the three strains in 2D by 



(la) 
(lb) 
(lc) 



ei = (?/n + r?22)/v2 , 
e2 = (mi - te)/v2 , 
e3 = {Vi2 + ??2i)/2 , 

where the components of the strain tensor 77 are 
1 / 

Vij = 2 + U it + u kdUk,j) 5 



(2) 



here u 



1,3 



dui I dxj 



and repeated indices are 



summed. All three strains vanish in the T phase. The 
deviatoric strain e2 is the primary order parameter of the 
T-0 transformation. In the lowest-energy product state, 
e2 takes one of two degenerate values ±e2o corresponding 
to a stretch in either the X\ or the x 2 direction. The 
dilatational and shear strains e\ and ^3 vanish for these 
two states, and also for twin bandsEl, but not for the 
complex domain patterns formed by colliding bands. 

At this stage in the theory of ferroelastics, one wants to 
examine the simplest possible form for the energy density 
T , to include only those terms required by symmetry, for 
stability, and to explain experiment. We start from the 
expression 



T 



A 



A 2 (T)e 2 2 



A 21 , B 2 4 C2 6 



+ y (Ve 2 ) 



(3) 



all terms are invariant under the symmetry operations 
of the T group. The dilatational, deviatoric and shear 
stiffnesses Ax, A 2 (T) and A3 in the first term are re- 
lated to the elastic constants. Stability requires A\ > 
and ^3 > 0. But A 2 (T) softens with decreasing T, as 
A 2 (T) = a(T - T ), and the T phase is unstable for 
T < Tq. To describe the phase transition, we need the 
terms in e\ and e\\ we assume a first-order transition 
(B 2 < 0), and so C 2 > for stability. At high T, namely 
A 2 (T) > B\j\C 2 , only the T minimum exists. At lower 
T, two O minima occur at e 2 = ±e 2 o(T), where 



e 20 (T) 



-B 2 + 



^Bi-4A 2 (T)C 2 ) /(2C 2 ) 



1/2 



(4) 



At the transition temperature T c , found from A 2 (T C ) = 
3S|/16C2, the three minima e 2 — 0, ±e 2 o(T c ) are degen- 
erate; here e 2 o(T c ) = \J~-iB 2 jAC 2 . Finally, the gradient 



term is responsible for the wall en^gic-the other deriva- 

largely because 



tive invariants^ are unimportant ,rfll§iEf ' 
the primary physical spatial dependence is in e 2 . 

The parameters of the theory are not well known for 
any material. To reduce the number of unknown param- 
eters, and possibly obtain a universal theory that applies 
qualitatively to many materials, we transform variables 
by 



Bj -» [e 20 (T c ) x 10 3 ] e 3 , 
%i — * \J d 2 /A 2 (T c ) Xi , 
T^A 2 {T C ) [e 20 (T c ) x 10 3 ] 2 T ■ 



(5) 
(6) 
(7) 



also, we define the dimensionless temperature r = 
A 2 (T)/A 2 (T C ) — (T — T )/(T C - T ) and dimension- 
less stiffness parameters £i = Ai/A 2 (T C ) and C3 = 
A 3 /A 2 (T C ). The scale factor in Equation (||) is chosen 
so that the deviatoric strain at T c is 10 -3 , an arbitrary 
value; the hidden but necessary assumption here is that 
the strains are small and so the nonlinear term in Equa- 
tion (^J) can be neglected. The energy density in terms 
of the new variables is 



1 



6 



^(Ve 2 ) 2 (8) 



where b = —4 x 10 6 and c = 3 x 10 12 . If temperatures 
near T c are accessible, the three parameters in Equation 
(^) can be determined from the elastic constants just 
above T c , the strain e 2 at T c and the T-dependence of 
e 2 . For YBa 2 Cu3C>7, typical values at low T arec3 an or- 
thorhombic distortion of 2(6 — a) /(a + b) = 0.017 (giving 
620 = 0.012), and a wall width of « 1.3 nm. 

Static structures predicted by Equations (||) (or (g)) 
are discussed in Refs. |(|||,[l2| Domain walls have lowest 
energy (ei and e% are zero) when in the T 110 and 110 
planes. The walls link the variants but also rotate them 
by an angle proportional to e 2 . The rotation, which has 
no counterpart in conventional order-parameter systems, 
gives rise to unusual effects when orthogonal walls collide; 
for example, the visual wall length increases in the col- 
lision region, due to variant narrowingE3 resulting from 
formation of a disclination. 

Different structures are found in soft or stiff systems, 
depending on whether the energy cost is small or large for 
wall directions off the optimal 110 and 110 planes. The 
relevant parameters are the ratios C1/C2 and C3/C2 of the 
dilatational and shear stiffnesses to the deviatoric stiff- 
ness £2 = T+36e2o+5ce2 . The energy cost increases with 
both ratios, though more strongly with (1/C2 it seems. 
Strangely, systems soften with decreasing T, because £2 
increases ((2 — ► 4 as r — ► 1 from below, with £2 = 238 
at r = — 50 for example). If a system is moderately stiff 
just below T c , then features like split tips characteristic 
of stiff systems may disappear on cooling, provided that 
low enough temperatures are accessible and that the re- 
laxation is not too sluggish. 
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B. Time evolution 



The Lagrangian density is 



1 



C = T-V= -p(u 



T 



(9) 



where T is the strain-energy density. To represent the 
nonconservative forces in the system, we use a Rayleigh 
dissipative functional with density 



n 



1 



(10) 



here ij = dej/dt. This form respects the symmetry of 
the T phase; it assumes evolution without plastic flow. 
The important point is that Equation ( |ici| ) leads to dissi- 
pative forces that are functions of the spatial derivatives 
of the velocity, as one would expect, since uniform mo- 
tion of the material canaot dissipate energy. Then the 
equations of motion are£3 



pUi 



ik,k 



&ik.k 







where 



a ki 



dK 

dui,k 
dT 
du lyk 



(11) 

(12) 
(13) 



We assume that the dissipative term is much larger 
than the inertial term. This approximation fails however 
at small wavenumber, as discussed for example in Ref. 
|l3| . In particular, the zeroth Fourier component should 
be considered separately since the last two terms of Equa- 
tion ([ll]) are then zero; then the inertial term tells us that 
the motion is uniform, determined by the initial value. 
Without the inertial term, Equation (111]) simplifies to 



ik,k 



'Oik.k 



(14) 



The summation on the index k prevents integration of 
these equations, except in ID. In ID, the constant of inte- 
gration is crucial, for it represents external forces applied 
to the boundary that may hold the system in a static 
configuration that is not necessarily the unconstrained 
minimum of the strain energy. 

The equations of motion ( |T4| ) in terms of the strains 
are 



A[e hl + A' 2 e 2 ,i 
A\t\, 2 - A' 2 e 2 ,2 



e3,2 



-(Gi 



G 2 



A' 



G 



:i:2) 



G 3 ,i) 



(15a) 
(15b) 



where Gi — ST/Sei and the individual functionals are 

Gi = A x e x , (16a) 
G 2 = A 2 e 2 + B 2 e\ + C 2 e\ ~ d 2 V 2 e 2 , (16b) 

(16c) 



G 3 - A 3 e 3 /V2 



We emphasize that our equations of motion (15) are 
not those of time-dependent Ginzburg-Landau (TDGL) 
theory. Schematically, the latter are 



-bT /Set 



(17) 



with a nonlocal expression^ for the density T\ the ma- 
jor difference is the additional space derivatives on both 
sides of Equation (15). Equation ( jl7| ) has much intu- 
itive appeal, not least because it continues the analogy 
with fcrromagnets. Nevertheless, it cannot be correct in 
principle, and in fact its predictions disagree with those 
of Equation (15). We illustrate the point by consider- 
ing a material with short-range internal forces, uniformly 
stretched by external forces applied at the ends. When 
the latter are abruptly released, relaxation begins at the 
ends and propagates inward, taking a finite amount of 
time to reach any point in the bulk; the ions (except 
those near the ends) feel equal but opposite forces from 
their neighbours until the disturbances reach their vicin- 
ity. Equations (15) have the correct behaviour, whereas 
Equation (|l7]) predicts instantaneous response. 
Equations (15) differ also from the dynamics 



iii cx —ST/5ui 



(18) 



of Refs. ^ and |T^, the former at T — 0. Not having exam- 
ined physical settings comparable to those where Equa- 
tion ( Jlq ) was used, we cannot compare its results with 
those of Equation (15). The right-hand side of Equation 
( p8| ) agrees that of Equation (jl4|); but the left-hand side, 
a dissipative force proportional to the velocity, cannot be 
correct in principle. 

From Equations (15), the equations of motion for the 
two components of u are 



A' 2 )dl + m/2 (A[ 
A' 2 + 4/2)did 2 (A[ + A' 2 )d 2 2 + A' 3 dl/2 



A' 2 + ^72)5^2 



G\.\ 
G\ 2 



G 2t i 
G 22 



G 3 , 2 
G 3 .i 



Ml 
»2 

(19) 



By solving these equations, we satisfy automatically the 
2D compatibility relation 



V 2 ei 



Of 



d 2 2 )e 2 



V89i5 2 e 3 - 



(20) 



in the small-strain approximation. This necessary and 
sufficient requirement that the strains be derivable from 
the displacement can be obtained by starting from u^i2 = 

M»,21- 

The three viscosity parameters A\ are not known from 
experiment, though of course all must be > 0; it is then 
reasonable to consider the simplest possible theory. In 
choosing parameter sets, we should avoid those that give 
a vanishing determinant 



1 



Det= -(A' 1+ A' 2 )A' 3 (dt 



a 9 4 ) 



[^[A'.-iA'.-A'^dfdi 



(21) 
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of the coefficients on the left-hand side of Equation (19) 
inspection shows that only one of the A\ can vanish 
Other cases of interest are those for which the determi- 
nant factors, i.e. 



AA[A' 2 - (A[ - A' 2 )A' 3 = ±(A[ + A' 2 )A> 3 



3 ) 



(22) 

giving three possibilities: (a) A[ = 0, (b) A 3 = 2A' 2 , and 
(c) A' 2 = 0; a fourth, namely A 3 = —2A[ fails on grounds 
discussed above. Since e2 is the primary order parameter, 
we should keep A' 2 ; the time scale is then adjusted so that 
A' 2 = 1. The choices A' 3 — 2A' 2 (the isotropic case) and 
A\ = are convenient, for then the left-hand sides of 
Equation ( |l9| ) decouple. We verified that taking A[ = 1 
vs. A[ — has little effect during the evolution; the fully 
relaxed configurations can differ however. 

We imposed periodic boundary conditions on the dis- 
placement u, thereby forcing domain walls into the sys- 
tems; the equilibrium states are a single twin band, opti- 
mally with only a pair of walls. We solved Equations ( |l9| ) 
using a finite-difference, fast-Fourier-transform method. 
At the beginning of each time step, the displacement 
field was known at each point of the space grid. Finite- 
difference approximations (centered on a 5 x 5 grid) were 
used to compute the derivatives and so to obtain the 
right-hand sides in real space. The latter were then 
Fourier transformed. The Fourier components on the 
left-hand sides were found using the same finite-difference 
approximations and then advanced in time using the Eu- 
ler method (with time step 10 -5 or so). The results were 
then Fourier transformed back to real space to begin the 
next step. 



III. T-O NUCLEUS IN TWO DIMENSIONS 

This section studies the nucleus resulting from perturb- 
ing the supercooled T phase in various ways. All results 
are for a grid of 512 x 512 points, with step size 0.4. 

We first present results obtained by displacing a sin- 
gle point off a high-symmetry direction. Figure 1 shows 
snapshots for r = — 50 (C2 = 238) and for four sets of 
values of £i and £3, all at time t = 0.18 after identical 
nucleation events; the viscosity parameters are A' x = 0, 
A' 2 = 1 and A 3 — 2. Figure 2 shows snapshots of the 
same systems at the later time t — 0.24. Very little is 
known about the relative importance of the stiffnesses £1 
and £3 and so we investigated some extreme cases; we 
find stronger dependence on £1 than on £3. Parts (a) and 
(b) of Figures 1 and 2 show soft systems (Ci = 1), with 
£3 = 1 and 1000 respectively, whereas parts (c) and (d) 
show moderately stiff systems (£i = 1000), again with 
£3 = 1 and 1000 respectively. The important point is 
that the nucleus has very different shapes in soft and 
stiff systems; one notes also the more rapid growth in 
the latter. 

In the soft systems, the domain walls lie off the opti- 
mal directions; the nucleus retains its disk shape as it 
expands. 



In the stiff systems, the domain walls are much closer to 
the optimal orientations. The nucleus has a striking X 
shape with arms in the 110 and 110 directions; growth 
transverse to the arms results from the appearance of new 
variants near the nucleation site and their subsequent 
growth along the arms. 

Other sets of simulations started from point displace- 
ments in high-symmetry directions (100 and 110), and 
others from displacements of small areas. Every soft sys- 
tem gave a disk with eight or more distinct domains em- 
anating from the disturbed area; every stiff system gave 
the X-shaped nucleus. 

Yet more sets of simulations started with disk-shaped 
regions containing several parallel stripes in one direc- 
tion, this in a bid to approximate the nucleus reported 
in Ref. [H]. These attempts gave nuclei much like those 
from a point perturbation. Parts (al) and (a2) of Figure 
3 show that a soft system evolves toward the flower-like 
patterns in parts (a) and (b) of Figures 1 and 2. Parts 
(bl) and (b2) of Figure 3 show the evolution of a stiff 
system; the overall size of the figures is identical to those 
for the soft system, but the area of the starting configura- 
tion is about one- fourth that in parts (al) and (a2). Fast 
growth occurs parallel to the starting walls, but twinned 
jets shoot out in the transverse direction, thereby evolv- 
ing the system toward the X shape in parts (c) and (d) of 
Figures 1 and 2. The two sets of jets are more asymmet- 
ric here, because the rapid longitudinal growth exagger- 
ates the greater asymmetry in the starting configuration. 
Nevertheless, it is clear that even this starting configura- 
tion is also unstable toward the formation of perpendic- 
ular jets and evolution to the X shape. 

Simulations at other temperatures (between r = —100 
and r = —5) gave results qualitatively similar to those 
described in Figures 1 to 3; the major difference is that 
the nucleus grows more slowly at higher T, as expected. 
The important point is that the flower/X shapes were 
found for soft/stiff systems at all T. We were unable to 
nucleate the low-T phase above r — — 5 (well below the 
stability limit r = of the T phase) and so we could not 
examine the parameter set of Ref. |t| 

Because the gross features are independent of the start- 
ing configurations and temperature, we believe that we 
have found the nucleus of the T-0 transformation in 2D, 
with possible application to thin films. It is reasonable 
to expect that x-y cuts through the 3D T-0 nucleus will 
resemble our 2D nucleus. 

None of our simulations (with any starting configura- 
tion, with either soft or stiff parameters, at any tempera- 
ture) gave a nucleus resembling that found using TDGL 
theory in the strains. The 2D T-0 nucleus of Ref. [ll] ac- 
cords with one's intuition based on conventional systems. 
It is compact, elliptical in shape (with axes along the 110 
and 110 directions), and internally twinned (with walls 
parallel to the major axis); the twinning generates both 
positive and negative displacements which largely cancel 
overall. Transverse growth occurs by adding walls and 
variants, whereas existing variants grow only longitudi- 
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nally. Although other aspects are different (Ref. |Tl| stud- 
ied soft systems, used a somewhat different strain-energy 
functional, and worked at higher T, namely r = 0.3), 
it is likely that the different results reflect the different 
dynamics. 

None of our simulations gave a nucleus like that in the 
more phenomenological study of Ref. [2^, namely growth 
to an untwinned square which then flowers. 

The only previous use of the equations of motion 
( |l4| ) to exapine nucleation was in a study of H-0 
ferroelastics;tZl these systems are dominated by disclina- 
tions. In soft systems the nucleus is flower-like, as in T- 
O systems, but has 12 arms; in stiff systems it branches 
early in the growth, without forming the long arms seen 
above in T-0 systems. 



IV. COARSENING 

This section studies the coarsening phenomena that oc- 
cur after completion of the phase transition. The interest 
lies in the unconventional behaviour relative to that ob- 
served in order-parameter systems. Simulations started 
from systems with orthogonal twin bands, relaxed inter- 
nally but not in the collision regions. The initial relax- 
ation from these artificial high-energy configurations is 
rapid and of no interest; we present results at later times, 
but well before equilibrium is reached. 

Figure 4 shows four pairs of snapshots. 
Parts (a) to (c) are for soft systems with different initial 
conditions, all with parameters A[ = A' 2 = A' 3 = 1, 
Ci ~ C'3 — 10, and t = —50; the times between the pairs 
are 0.5, 0.5 and 1.0 respectively. 

In part (a), the island at the centre vanishes, but other 
islands form as some narrow domains pinch off and re- 
tract. 

In part (b), one tip retracts to form rank with its neigh- 
bour; at the lower right, other tips retract in unison, 
keeping the rank. 

In part (c) , coarsening occurs by different kinds of coordi- 
nated events; domain merges parallel to the smaller-scale 
patterns occur at the top left and perpendicular at the 
bottom right. 

Part (d) corresponds to a stiffer system, with parameters 



A' 



1, A' 



2, Ci = Cs 



500, and r 



-100 



(£2 = 452); the time difference is 0.6. The patterns are 
strikingly similar to Figures 7.9 and 7.17(b) of Ref. | 
and to a lesser extent Figure 2(b) of Ref. 3(b). One 
sees the formation of a split tip and also the counterin- 
tuitive variankjiarrowing and wall wobbling found in the 
static theoryJl3 Related theories of needle twins and tip 
splitting are given in Rcfs. 
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The observation of tip splittingd'cl in YBa2Cri307 sug- 
gests that this material is moderately stiff (Ci ^ C2) a t 
the temperatures investigated. Values of the elastic con- 
stants suggest that Pg-Pd alloys (cubic-tetragonal) are 
also moderately stiff.E3 



These coarsening phenomena, like the nucleation phe- 
nomena reported in Section 3, confound intuition based 
on conventional order-parameter systems. The relax- 
ation cannot be characterised by any simple rules; that is, 
the changes from one snapshot to the next cannot be pre- 
dicted by inspection of the strain patterns alone. The vis- 
ible domain^wall length often increases. The relaxation 
is nonlocal;E!j rapid changes occur in one part of the sys- 
tem while other parts, with no apparent major differences 
from the first, stay almost unchanged. The tendency is 
toward coarser patterns, but occasionally the topology 
becomes more complicated (as when islands form). The 
ribbons seldom retract immediately, even though retrac- 
tion reduces the wall length. Particularly odd are the 
rank formation of tips and their linked withdrawal, the 
variant narrowing and the splitting of tips. Transverse 
wall motion occurs only locally, for example in the pro- 
cess of pinching off the other variant. 

Our simulations resemble in some respects those of 
Refs. HH and less those of Refs. [j|]||ll]. 

Coarsening mechanisms in simulations of H-0 



sy stems, O also using Equations (14), differ from those 



in Figure 4 (again due to the disclinations in H-0 sys- 
tems). 



V. SUMMARY 

We have derived general equations of motion for proper 
T-0 ferroelastics including inertia, dissipation and inter- 
nal elastic stress. These equations, and more importantly 
their predictions, differ from those of all previous stud- 
ies of proper T-0 ferroelastics. We studied the growth 
of the O nucleus for both soft and stiff systems, in 2D. 
The soft system expands as a disk with time, while the 
stiff system assumes a characteristic X shape, with twin- 
ning along the arms. We studied also the coarsening 
mechanisms that relax the O phase toward local equilib- 
rium, again in 2D. We observed the formation and dis- 
appearance of island domains, tip retraction and domain 
merging, both parallel and perpendicular to existing do- 
main walls; in stiff systems we observed the formation 
of split tips. These mechanisms are likely not observ- 
able in proper ferroelastics, because the time scale is ex- 
pected to be short; likely one can examine only patterns 
in quenched samples. Perhaps they are observable in 
improper systems, where the time scale may be longer; 
again, our strain-only theory does not apply in principle 
to improper ferroelastics, but it explains many puzzling 
features of patterns reported in Refs. || and ||, and so 
perhaps it can shed light on the dynamics also. 

The above treatment should be extended to include 
thermal noise, first to examine the early stages of nucleus 
formation, and second to allow the system to surmount 
energy barriers. The tweed structure should be examined 
in the presence of noise, perhaps also with compositional 
fluctuations. The incrtial term should be examined to 
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determine whether it affects the dynamics significantly. 
The difficulty is to find realistic values of the viscosity 
parameters; one can easily be misled here. 

The primary need in the field is however in situ obser- 
vations of the dynamics in T-0 systems; these are diffi-. 
cult and correspondingly rare. The available studiesErlEl 
cannot decide the relative merits of the many theories. 
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FIG. 1. Greyscale snapshots of orthorhombic (O) nu- 
clei growing after identical perturbations of the supercooled 
tetragonal (T) phase. The two O variants are white and black, 
the T matrix grey. The four parts correspond to different 
choices for the dilatational and shear stiffnesses; the parame- 
ter values are given in the text. 
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FIG. 2. The nuclei of Figure 1 at a later time. 




FIG. 3. Snapshots of O nuclei growing in a T matrix. The 
starting configuration was a disk, internally twinned. Parts 
(al) and (a2) show the nucleus at times 0.12 and 0.24 for a 
soft system (A[ = 0, A' 2 = 1, A' 3 = 2, £i = Cs = 1, r = -50). 
Parts (bl) and (b2) show the nucleus at times 0.12 and 0.18 
for a stiff system (A[ = 0, A' 2 = 1, A' 3 = 2, £i = (3 = 1000, 
r = -50). 
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FIG. 4. Pairs of snapshots showing the time evolution of 
structures for four different initial conditions or parameter 
sets. Each part is a 128 x 128 piece of a full 256 x 256 sim- 
ulation with step size 0.2. Parts (a) to (c) are soft systems 
and part (d) stiff; the parameter values and time intervals are 
given in the text. 
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